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ABSTRACT. A parameter-dependent model involving nonlinear diffusion for an age-structured population is 
studied. The parameter measures the intensity of the mortality. A bifurcation approach is used to establish 
existence of positive equilibrium solutions. 



1. Introduction 

Let u = u(t, a, x) > be the distribution density at time t > of individuals of a population structured 
by age a G J := [0, a m ] and spatial position x G f2, where a m G (0, oo) denotes the maximal age and £1 is 
a bounded and smooth domain in Mr. Suppose that the individual's movement is governed by a nonlinear 
diffusion term div x (D(U(t, x), a)W x u) with dispersal speed D(U, a) > depending on age and on the 
total local population 



U(t,x):= / u(t, a,x)da. 
Jo 

Let b = b(U, a) > and ft — p,(U, a) > denote respectively the density dependent birth and death rate. 
Then a simple model describing the evolution of the population with initial distribution u° = u°(a, x) > 
is 

dtu + d a u — dxv x (D(U (t, x), a)V x uj — jl(U(t, x), a)u , i>0, a£ J, xgO, 

u(t, 0, x) = I b(U (t, x), aju(a) da , t > , X G Q , 
Jo 

5u(t, a, x) + (1 — 5)d v u(t, a, x) = , OO, b£ J. ig dfl , 

u(0, a, x) — v°(a, x) , a G J , x G £1 , 

where either S — 1 or 6 — corresponding to Dirichlet or Neumann boundary conditions. Models of 
this type have a long history and we refer to ll22l for a survey of structured population models. The well- 
posedness of these equations and related population models involving nonlinear diffusion is investigated 
e.g. in Ifl9ll . Questions regarding the large time behavior are linked to equilibrium solutions. In this paper 
we thus shall focus on nontrivial nonnegative equilibrium solutions for such equations, that is, on time- 
independent solutions u = u(a, x) > with u ^ 0. 

Positive equilibrium solutions for age-structured population models without diffusion are studied e.g. in 
(9] [TOj [TTj using bifurcation techniques or also in lfT31l using fixed point theorems in conical shells. A bi- 
furcation approach to age-structured population models with linear diffusion and linear birth but nonlinear 
death rates is used in |fT3l[T4ll . For an approach to age-structured models including both nonlinear diffusion 
and nonlinear death and birth rates we refer to ||20ll2Tl . where local and global bifurcation, respectively, is 
shown for a bifurcation parameter measuring the intensity of the fertility similarly as in [9, 10, 1 1 ]. The aim 
of this paper is to demonstrate that also the intensity of the mortality can be treated as bifurcation parameter 
in age-structured models with nonlinear diffusion. Moreover, as expected and opposed to the results of e.g. 
where the fertility intensity varies, in the present situation subcritical bifurcation occurs under realistic 
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assumptions. 

The approach we choose is based on introducing a parameter A measuring the intensity of the mor- 
tality without changing its structure; that is, we shall consider parameter-dependent death rates of the 
form p, — \fi(U, a) with /i = fJ,(U, a) being a fixed reference function. Thus we are looking for solutions 
u = u(a, x) to the parameter-dependent problem 

d a u = div x (D(U(x),a)V x u)—Xfi(U(x) 7 a)u, a£ J, i£fi, (1.1) 



u(0,x)= b(U(x),a)u(a)da , x E Q , (1.2) 

Jo 

Su(a, x) + (1 — 5)dvii(a, x) = , a e J , x 6 dfl . (1.3) 

Clearly, u = is a solution to dl.lt - dOT l for any value of A. The main goal is then to establish existence 
of nontrivial solutions which are also nonnegative. Under suitable assumptions we shall prove that the 
theorem of Crandall-Rabinowitz 1 8 1 applies so that there is a unique value Ao > for which a nontrivial 
branch {(A, u\); | A — Ao | small} bifurcates from the trivial branch {(A,0); A G R} at the critical point 
(Ao,0) and that at least one part of the nontrivial branch near the critical point consists of nonnegative 
solutions. 

To be more precise, let o\ be the first eigenvalue of —A x on subject to Dirichlet (if 6 = 1) or Neumann 
(if (5 = 0) boundary conditions, hence o\ > in the first and o\ — in the second case. Suppose that 

/ 6(0,a)e~ CTl -'° aD(0 ' r)dr da > 1 and /x(0, a) > for a near . (1.4) 
Jo 

Roughly speaking, the first assumption in dl.4t may be interpreted as that for a zero death rate, the popula- 
tion is (locally) increasing. Letting Ao > be such that 

' 6(0, a) e- Xo f° ^°' r)dr e- ai f° D (°' r ) dr da = 1 , 

the following result on local bifurcation holds for equations dl.lt - dl.3t : 

Theorem 1.1. Let D e C 00 - 1 ^ x J) with D(z,a) > d > for z G R and a G J. Further, let 
/i, b G C 00,1 (R x J) be nonnegative and suppose ( 1 1 ,4b . Then (Ao, 0) is a bifurcation point for dl.lt - dl.3t , 
that is, there is a unique local branch of nontrivial nonnegative solutions 

(A,u) in r + x (c(j,c{n)) nc^j.c^)) nc(j,c 2 (n))) 

emanating from the critical point (Ao, 0), where J :— J \ {0}. In addition, if 5 — and 

b(z, a) < 6(0, a) , fj,(z, a) > p,(0, a) , z > , a G J , (1-5) 
f/ien bifurcation is subcritical, i.e. A < Xofor any nonnegative solution (A, w). 

Assumption dl.5t is a common modeling assumption stating that effects of population densities do nei- 
ther increase fertility nor decrease mortality. The result thus shows that lowering the intensity of mortality 
below a critical value leads to other equilibrium solutions than the trivial one. We also refer to Section[3]for 
an example where subcritical bifurcation occurs when 6=1 and dl.5t holds. 

We shall emphasize that Theorem ll.ll is merely a consequence of the considerably more general Theo- 
rem l2.5l that includes general nonlinear elliptic diffusion operators not necessarily in divergence form (and 
also less regular data). The proof of Theorem 1 1.1 1 will be given as an application of Theorem 12. 5 1 in Sec- 
tion [3] 

To cover a great variety of applications we thus shall consider dl.lt - dOl l as an abstract equation of the form 

d a u + A(u, a) u = —\h(a)u + g(X, u, a)u , a G J , (1-6) 

u(0) = / b(u,a)u(a)da, (1.7) 
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in an ordered Banach space Eq with positive cone Eq for the unknown function u : J — > Eq. Here, 
A(u, a) defines for fixed (u, a) G Eo x J a bounded linear operator from a subspace E\ of Eq into Eg. 
Problem (ll.U - dl.3b then fits into this abstract framework by choosing 

£0 := L q (Q) , Ei := {v e W^(Q); 5v + (1 - S)d v v = on dQ} 

for some q G (1, oo) (where boundary values are interpreted in the sense of traces) and letting 

A(u, a)w := — div x (D(U, a)V x w) , fo(a) := /i(0, a) , and <?(A, u, a) := \(n(U, a) — /z(0, a)) . 

In Section|2]we consider the abstract equations dl.6l ). ( 11.7b and prove under suitable assumptions in Theo- 
rem |2.5| a local bifurcation result. In Section 3 we give applications of Theorem |2.5| and prove in particular 
Theorem 11.11 Finally, the appended Section [4] contains a result on the differentiability of superposition 
operators in Sobolev-Slobodeckii spaces used for the applications in Section[3]that we were unable to find 
in the literature in this form. 

2. The abstract bifurcation result 

Studying the nonlinear problem (11.61 . ( 11.7b demands an investigation of its linearization around u = 0. 
We first state the precise assumptions required. 

2.1. Assumptions. Given Banach spaces E and F we let C(E,F) denote the space of all bounded and 
linear operators from E into F, and C(E) :— C(E, E). We write Cis(E, F) for the subspace of C(E, F) 
consisting of all topological isomorphism and fC(E, F) for the subspace of compact operators. 

For the remainder of this section let J := (0, a m ) with a m G (0, oo] and note that J may be bounded 
or unbounded. Moreover, we fix an ordered Banach space Eq with positive cone Eq and a dense subspace 
Ei thereof which is also supposed to be compactly embedded in Eq. This latter property we express 
by writing Ei s> Eq. Given 9 G [0,1] and an admissible interpolation functor (•, -)g we equip the 
interpolation space Eg := (Eq, Ei)g with the order induced by the positive cone E^ := Eg n Eq . Note 
that E-0 Eg for < 9 < $ < 1 according to J5] I.Thm.2.11.1]. In particular, we fix p G (1, oo) and 
set E^ := (Eo, Ei)i_i/ p p ; that is, E q is the real interpolation space between Eq and Ei of exponent 
<; := 1 — 1/p. We then assume that 

mt(E+) ? , (2.1) 
where mt(E^) denotes the topological interior of the cone E+, We set 

E := L p (J, E ) and Ej := L p (J, E x ) n W$(J, Eq) 

and recall that Ei BUC(J, E q ) (see 0). Thus, the trace operator 70U := u(0) for u G Ei is a 
well-defined operator 70 G £(Ei,i? f ). We also set := Ei n L+(J, Eq). Suppose that 

F is a Banach space ordered by a positive cone F + with F ■ E\ ^ E f and F + • Eq ^ Eq , (2.2) 

where e.g. F ■ E\ -E ? means a continuous bilinear mapping (i.e. a multiplication) F x £1 — > E 1 ^, 
(/, e) 1— ► / • e. Let E be a fixed ball in Ei centered at of some positive radius and assume that 

g G C 1 (M + x E, L^J, F)) with g(X, 0) = for A G K+ (2.3) 

and 

h G L+( J, R) n Loo (J, R) with /i > near a = . (2.4) 

Observe that ( 12.3b guarantees that we may interpret g(X, u) as an element of L^J, C(Ei, Eq)) for (A, u) 
in IR + x E fixed. Suppose then that 

AeC\i:,L 00 (J,C(Ei,EQ))) (2.5) 
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is such that 

A(u) + Xh — g(X, u) £ L^J, £{Ei,Eq)) generates a positive parabolic 

evolution operator II(a,u) (a, cr), < a < a < a m , on Eq with regularity subspace (2.6) 

£a for each (A, ti)el + xE. 

We refer to [5] for a definition and properties of parabolic evolution operators. Note that, due to ( 12. 3t and 
(12.41 , the parabolic evolution operator Ilo := II( ) is simply generated by ^4(0) and 

n (A ,o)(o, a) = e- x ^ h ^ dr IL (a, a), < a < a < a m . (2.7) 
We further suppose that 

there are £ E R, p, w > 0, k > 1 such that C + A(0) E C p [ {J,U{E 1 , E ; k,u)) (2.8) 

and that 

A(0) + Xh £ ioo(^, £(Ei,Eo)) possesses maximal L„ -regularity on J , 

(2 9) 

that is, (d a + A(0) + A/1,70) € £«s(Ei,E x£ t ), for each A > . 

We refer again to [5| for a definition of the space %(Ei, Eq; k,uj) and details about operators having 
maximal L p -regularity. We agree upon the notation A(u, a) := A(u)(a) and e.g. (hu)(a) := h(a)u(a) for 
a e J and u E Ei. We point out that, owing to ( 12.61 ). ( 12.91 ), and [5, IILProp. 1.3.1], the linear problem 

d a u + (A(0, a) + Xh(a)) u = f(a) , aeJ, u(0) = u° 

admits for each datum (/, u°) E Eo x E q and A > a unique solution u E Ei given by 

u(a) = n (Aj0) (o,0)t»°+ / n (A , 0) (a,(7)/((7)ckr, aEJ, (2.10) 
Jo 

satisfying for some c = c (A) > 

\\u\\ El <co(||/||E + ||u b t ) ■ (2.11) 

Moreover, invoking [5, II.Lem.5.1.3] it follows from ( 12.8b that there are Mq > 1 and ljq E R such that for 

0<7<^<a<l 

l|n (a,a)|| £(£7) + (a- ( 7) Q ^||n ( a , ( 7)|| £( ^ iB£>) <M oe -«( a - CT ) , < a < a < a m . (2.12) 
We also assume that 

n (a, 0) is strongly positive for each a E (0, a m ) , (2.13) 

that is, n (a, 0)0 E Lat(J5+) for <j> e E+ \ {0} and a £ (0, a m ), and that 

II (a, 0)n (cr, 0) = II {a, 0)n (a, 0) , < a, g < a m . (2.14) 

The latter condition means that the operators {A(0, a); a £ J} commute with each other. Finally, we 
assume that 

b e C l {V,L+{J,F)) withO ^ b Q := b(0, •) £ L+(J,R) 

/■<^ ' (2.15) 

and / b (a)e u ° a da < 00 , 
Jo 

where p' is the dual exponent of p, i.e. 1/p + l /p' = 1. The last condition in (12.15b is obviously superfluous 

if a m < 00 or luq < 0. 

For the remainder of this section we assume that conditions (12. lb -( f276T >, (12. 8t , (12. 9K (12.13l l-( |2~T5T l hold 
and refer to Section[3]for examples where these conditions are met. In particular, they hold in case of The- 
orem ll.ll 
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2.2. The linear problem. We begin by investigating the linearization of ( 11.61 ), ( 11.7b around u = 0, that is, 
by investigating the problem 

d a u + (A(0,a) + Xh(a))u = 0, a e J , (2.16) 

u(0) = / 6 (a)"(a)da . (2.17) 
Jo 

It readily follows from the previous observations that any solution (A, u) G K + x Ej of ( I2.161 l, ( 12.17b is of 
the form 

u (a) = e- A ^' l W d '-no(a,0)u(0) , a e J , u(0) = Q A u(0) , (2.18) 
where the operator Qa is given by 

Q A := / Om 6o(o)e- x /o a ' l ( p > dp llo(o ) 0)da 
Jo 

and enjoys the following properties: 

Lemma 2.1. For any A > 0, Q A G is strongly positive. Hence, the spectral radius r(Q\) > of 

Q\ is a simple eigenvalue of Q\ and of its dual operator Q' x with eigenvector B\ G int(_B+) and strictly 
positive eigenfunctional B' x G E' q , respectively. Moreover, r(Q\) is the only eigenvalue with a positive 
eigenvector. 

Proof. Let 6 G (0, 1 /p) and set := 6»+<r > ?. Then we derive from fl2J2l) and d27T3T ) that Q A G , JS^) 
for A > and thus, since ^ J5 5 by [5, I.Thm.2.11.1], we have Q A G A > 0. Hence the 

assertion follows from the Krein-Rutman theorem (e.g., see lfl2l Thm.12.3]) and assumption ( 12.11 ) provided 
we can show that Q\ G IC(E^) is strongly positive. To fill this gap let /' be any nontrivial element of 
the dual space E[ of with (/', 4>) E , > for G E+. Let ip G E+ \ {0}. Then it follows from Q 
Prop.A.2.7,Prop.A.2.10]and(|rBthat (/', n (a, 0)<^) Es > Ofora G (0,a m ) since n (a, 0)^ G mt(E+) 
by d213), and thus 

(f,Q\<p)E t = r m feo(a)e- A ^ /l( '' )d ''(/',n (a,0)^}^da > 
Jo 

owing to J2.151 ). Hence Q\<p is an interior point of E+ again due to Q Prop. A. 2. 7, Prop. A. 2. 10] and 
assumption J2.ll ). This yields the strong positivity of Q\. □ 

We assume in the sequel that 

r(Q ) > 1 • (2.19) 
Observe that ( 12.18b implies that u(0) is (if nonzero) an eigenvector of Q\ to the eigenvalue 1. If u is non- 
negative, i.e. u G E^~, then necessarily u(0) G and so r(Q\) = 1 by the previous lemma. The next 
lemma shows that r(Q\) is strictly decreasing in A. Hence, if ( 12.191 ) does not hold, there is no admissible 
(i.e. positive) value of A for which the linearized problem ( 12.16b . ( 12.17b admits a nonnegative nontrivial 
solution. The interpretation of the operator Q\ is that it contains information about the spatial distribution 
of the expected number of newborns that a population produces when the birth and death processes are 
described by b(0, •) and \h = A/i(0, •), respectively, and spatial movement is governed by A(0, •). Hence, 
at equilibrium these processes yield exact replacement. Roughly speaking, assumption ( |2.19b may be inter- 
preted as that the population subject to birth processes and spatial dispersal increases locally if A = 0, that 
is, if no deaths occur. 

Under assumption ( |2.19b , the following lemma guarantees the existence of a unique value Ao > with 

r(Qx a ) = l- 

The following auxiliary result uses the ideas of fl4l : 

Lemma 2.2. The mapping [A i-4 r(Q\)] : [0, oo) — > (0, oo) is continuous, strictly decreasing, and 
liirix-xxj r(Q\) — 0. In particular, there is a unique Aq G (0, oo) with r(Q\ ) = 1. 
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Proof. Since &o > with &o ^ and h > with h > near a = 0, it readily follows from d2.13b 
analogously to the proof of Lemma I27TI that Q x — is strongly positive for £ > A > 0, that is, 

(Qx - Qd<t> e wt(E+) , 4>eE+\{0}, £>A>o. (2.20) 

Given A > 0, let B x G int(££~) and B' x G £' be the eigenvectors and strictly positive eigenfunctionals 
introduced in Lemma |2~T1 Then, for £ > A > 0, we deduce from (12.20b that 

r(Q x )(B' x ,B^)^ = (Q' X B' X ,B^ E ^ = (B' X ,Q X B ( ) E ^ > (B' X ,Q S B^ = r(Q^) (B' x , B ( ) E ^ , 

whence r(Q x ) > r(Q^) so that [A M> is strictly decreasing. Next, let A > (the case A = is 

analogous) and consider a sequence (Xj ) such that < Xj — > X. Given e > sufficiently smal we may 
assume that < A — e < Xj < X + e for all j 6 N. Note then that ( 12.4b implies 

Qx-e B x < e s ^ Q x B x = e*" 1 r(Q x ) B x 

with \\h\\i denoting the Li-normof h. Since B x G E+ we derive 

r(QA)e £||/l111 > r(Q A _ £ ) (2.21) 

from [ 12] Cor. 12.4] and d2.lt . Conversely, we have 

QxB x+£ < e^Q x+£ B x+£ = e^r(Q x+£ )B x+£ 

and thus, invoking again [12, Cor. 12.4], 

e elMl r(Q x+£ ) > r(Q x ) . (2.22) 
Therefore, combining d2.21b . ( 12.221 1 and recalling that r(Q x ) is strictly decreasing in A, we obtain 

e -NU r (Q A ) < r (Q x+£ ) < r(Q Xj ) < r(Q A _ £ ) < e e "' 1 " 1 r(Q x ) . 

Letting e — > implies lim^oo r(Q x . ) = r(Q A ), whence the continuity of the function A h-> r(Q A )- 
Finally, the assumption that h > near a = together with (12.12b and ( 12.15b easily entails that 

< r(Q x ) < HQaIUck,) ->• , A^oo, 
from which the assertion follows in view of ( 12.19b . □ 

2.3. The nonlinear problem. To investigate the nonlinear problem ( 11.6b . ( 11.7b we apply the theorem of 
Crandall-Rabinowitz [8|. Clearly, the solutions (X,u) — (Xq + t,u) of (11.6b . ( 11.7b are the zeros of the 
function 

, . _( d a u + A{u)u + (X a + t)hu - g(X +t,u)u \ 
1 ' U) — \ u(0) - J a '" b{u, a)u(a)da J ' 

Assumptions (1272b , dOi d2~4l >, dOi and d2~T3T > imply that 

F : (-A , oo) x £ -> E x E, with 0) = , t > -A . 

Moreover, it is easily seen that all partial derivatives Ft, F u , and F tu exist and are continuous and that, since 
g(X, 0) = 0, the Frechet derivatives at (t, u) — (0, 0) applied to ip G Ei are given by 

n /n M ( d a (fi + {A + X h)ip \ 



where Aq := A(0, •), and 



r„,(0.l)) T -= ( ^ ) . (2.24) 
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Recall that B Xo £ int(J5+) with ker(l - Q\ ) = span{£ Ao }. Then (12791 ) implies that n (Ao o) (•, 0)S Ao 
belongs to E+ . Moreover, for 

^o(w) : = / bo(a)u(a)da, 



(K f)(a):= U (Xo>0) (a,a)f(a)da, f e E , 
Jo 

we deduce from d272"t , ([2TT0J, (|2TTj, and J27T3T ) that 

4 6£( E i>^). Kq&L^o^x). (2.25) 
With these notations we can state the following result. 

Lemma 2.3. L := F U (Q,Q) G £(Ei,Eo x is a Fredholm operator of index 0. In fact, 
ker(L) = span{n (Aoi0 )(-, 0)5 Aq } , 

rg(L) = {(<p, ^)6E xB f ;H € rg(l - Q Af) )} 

are both closed and dim(ker(L)) = codim(rg(X)) = 1. 

Proof. This is a reformulation of ll20l Lem.2. 1 ] using (I2.231 l. ( I2.251 l. (12.91 . Lemma lZTl and Lemma l2~2l □ 

This lemma also allows us to validate the transversality condition from (8). 
Lemma 2.4. We have F tu {0, 0)(n (Aoi0) (-, 0)B Xo ) £ rg(L). 
Proof. According to ( 12.241 ) and Lemma [23l we have to check that 

z e (K (hIl { x o , )(;0)B Xo )) # rg(l - Q Xo ) . 
Due to assumptions j2.4i , ( 12. 15b . and properties of evolution operators, we compute 



b (a) / n (AoiO )(a,CT)/i(cr)ri( Ao .o)(cr,0)B Ao dada 
Jo 

bo{a)(J h(a)do^j U {Xgfi) (a, 0)B Ao da . 

Thus z^O due to ( 12.41 ). ( 12.131 ). and (12.151 ). Using the commuting condition ( 12.141 ) we derive on interchang- 
ing the order of integration that 

Q Xo z = J 6 (s)n (ADi0 )(s,0) J b (a) (J h(a) do^j II (Aoi0) (a, 0)B Xo da ds 
bo (a) (^J h(a) do^j U (Xofi) (a, 0) J b (s) R( Xo ,o) (s, 0)B Xo ds da . 



Hence, simplifying the integral by recognizing Q Xo B Xa = B Xa in the integrand and reversing the compu- 
tations we obtain Q Xo z = z, that is, z £ ker(l — Q Ao ). But then z g" rg(l — Q\ ) since r(Q X(l ) = 1 is a 
simple eigenvalue of the compact operator Q Ao . □ 

Occurrence of local bifurcation in ( 11. 61 ), dl.7t is then a consequence of Lemma [2731 Lemma |274l and (8] 
Thm.1.7]. 

Theorem 2.5. Suppose d27lT)-<l276t. d278l> . d279l d27l3l - (l27l5t . d27T9b . Further let A > with r(Q Xo ) = 1. 
Then (Ao, 0) is a bifurcation point for ( 11.61 ), ( 11.7b . More precisely, there are Sq > and a unique branch 
{(A(e), u{e)) ; |e| < Sq} in R + x Ei emanating from (Ao, 0) with u(e) ^ if e ^ of the form 

u{e) = e(n (Ao , )(-,0) J B Ao + z(e)) , |e| < e . (2.26) 

Bof/z A : (— £o,eo) — > R + ana! z : (—So,£o) — ► Z are continuous, where Ei = ker(L) © Z wzf/z an 
arbitrary complement Z. Moreover, u(e) £ E^~ and Jqu(e) £ mt(E^) for £ £ (0,£o). 
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Proof. The existence of a nontrivial branch follows from Lemma 12.31 and Lemma 12.41 by applying (8] 
Thm.1.7]. It remains to prove the positivity assertion for e G (0,£o). Note that from (12.26b we have, 
for e G (0, £o) with £o sufficiently small, 

-7o«(e) = B\ +j z(e) G int(J5+) 

since z(e) — > in Ex ► BUC(J, 25 s ) and £?a G int(i/+). On the one hand, we derive 7ou(e) G and 
therefore, by positivity of the evolution operators assumed in (12.6b . 

= II(A( e ) itt ( e ))(-,Q)7ott(e) G E+ , £ G (0,£ ) . 

On the other hand, ^7om(£) and thus also 7ow(£) are quasi-interior points of E+, that is, (/, 7oit(£)} _e 5 > 
for each f e E'<-\ {0} with / > 0. So ( 12.11 ) and [7 Prop.A.2.10] imply that 7ow(e) is an interior point of 
E+. ' ' □ 

Remark 2.6. We proved a local bifurcation result for (11.6b , ( 11.7b under the assumption that h > 0. How- 
ever, the statement of Theorem \2.5\ still holds true if h < 0. 77ze only modification consists of replacing h 
by —h in assumption ( 12.4b so that the spectral radius r(Q\) is strictly increasing in A (see Lemma \2?2\ and 
one thus has, in addition, to replace ( 12.19b by the assumption that r(Qo) < 1. 

Moreover, ifh<0 one can even prove a global bifurcation result using the Rabinowitz alternative II 161 
provided that the nonlinearities in the operator A(u, a) are of "lower order", that is, if A(u, a) is a sum of 
operators ^4o( a ) + A*(u,a), where A (a) G £(Ei,E ) and ^4*(u, a) G C(Eg,E ) with 9 G [0, 1). The 
approach is similar to 121IL We also refer to [ 14] where the case h = — 1 is considered with linear diffusion. 

3. Examples 

Let f2 C M. N , N > 1, be a bounded and smooth domain lying locally on one side of dti. Let the 
boundary dil be the distinct union of two sets Tq and T% both of which are open and closed in dfl. Let the 
maximal age be finite, i.e. let a m G (0, oo) and set J := [0, a m ]. 

3.1. A general example. Consider a second order differential operator of the form 

A(U(x),a)w := -div x (D(U(x),a)V x w) +d(U(x),a) ■ V x w , (3.1) 
where, for some p > 0, 

D G C 5 ~' P (R x J) with D(z, a)>d o >0, z eR, a E J , (3.2) 

and 

d G C 4 ~' P (R x J, R N ) with d(0, •) = . (3.3) 
For simplicity we refrain from an explicit dependence of A on x G f2. Let 

^oGC^rx), m(a:)>o, leri, (3.4) 

and let v denote the outward unit normal to Y\. Let 



B(x)w 

Fix p, q G (1, oo) with 



to , on T , 

+ v Q (x)w , on Ti . 



2 A 

- + — < 1 , (3.5) 
P Q 



and let Eq :— L q :— L q (VL) be ordered by its positive cone of functions that are nonnegative almost 
everywhere. Observe that 

Ei ■= Wl B := Wl B {tt) := {u G W^ 2 ; Bu = 0} «-»• - E , 
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where Wi 2 (f2) is the usual Sobolev space of order 2 over L g (J7). Also note that, up to equivalent norms, the 
real interpolation spaces between Eq and E\ are subspaces of the Besov spaces B^ p := Bgf p (Q), that is, 

f B 2 q %, 0<2£<l/g, 
^ := (L„ W^ 2 B ) = ;8 := I {we B% ; U |r = 0} , l/q < 2£ < 1 + 1/g , 2£ 1 , 

[ {^e-B|| ;Bn = 0}, 1 + 1/g < 2^ < 2 , 

(see e.g. ifTTll ). In particular, due to ( 13.5b we have E s = B 2 q ~% v C 1+e (ft) for q = 1 - 1/p and some 
e > 0. So int(£+) ^ yielding <|270 , Fix any k G (2 - 2/p, 2) \ {1} and set F := ££ p;B with order 
induced by the cone of L q . Then pointwise multiplication 

is continuous according to ( 13.51 ) and [4, Thm.4.1]. Thus ( 12.21 ) holds. Let 

Ei := L P {J, Wis) n W^{J, L q ) and E := L P {J, L q ) 

and note that 

U := I u(a)da G Ex = W^ B , it G Ei . 
Jo 

Suppose that 

fi £ C 4 ~ :P (IR x J) , M>0, MO, a) > for a near . (3.6) 

Set h(a) := /«(0, a) and g(X, u)(a) := \(fi(U, a) - h(a)) for A G R, u G Ei, and a G J. Then ( 13.6b and 
Proposition l4.1l from the appendix ensure ( 12.31 ) and ( 12.41 ). Further suppose that 

beC 4 -' Q {Rx J) , b>0, 6(0,-) #0, (3.7) 

and set f>(u, a) := 6(C/, a) for u G Ei, a G J. Then ( 13.7b and Proposition l4.1l ensure ( 12.15b . Define 

A(u, a)u> := A(U, a)w , w G Ei , u G Ei . 

Proposition 14. 11 ( 13.2b , and ( 13.3b entail that the superposition operators induced by £>, d\D, and <i (again 
labeled D, d±D, and d) satisfy D, d±D, d G C^W^g, Loo(J, C 1+e (fi))). This yields 

whence ( [231 ). If u G Ei and A > are fixed, then A(u, •) + Xfi(U, ■) G C P (J, W(W£ B , L q )) from which 
we conclude ( 12.6b and ( 12.8b due to I.Cor.1.3.2, II.Cor.4.4.2] and the compactness of J. Noticing that 
A(0,a) = -D(0,a)A x by (E3 it follows from [l, Sect.7,Thm.ll.l] that for a G J and A > fixed, 
—A(0, a) — Xh(a) is resolvent positive, generates a contraction semigroup of negative type on each L r (f2), 
r G (l,oo), and is self-adjoint on £2(0). Hence A(0, •) + A/i possesses maximal L p -regularity on J 
according to J5] III.Ex.4.7.3,III.Thm.4.10.8], whence ( 12.9b . Moreover, since 

IT, (a, a) = ^(o.rJdrA. f < a < a < a m , 

where {e aAjc ; a > 0} is the semigroup associated with (— A x , B), condition ( 12.13b follows from the maxi- 
mum principle and (12.14b is obvious. Finally, let u\ be the first eigenvalue of (— A x , B) and let ip\ G W qB 
be a corresponding positive eigenfunction (e.g. see 0]). Then 

efo D(0,r)drA. ^ = e - ai /- D(0,r)dr v , x > < a < a m , 

and thus 

Qa<Pi = / b{ti,a)e- x f> {0 ' r)Ar D{0 ^ drA * Vl da = fc(A) ^i , 
./o 

where 

fc(A) := / 6(0, a) e~ x So ^°^ dr e^ 1 ^ ^r)Ar da _ 
Jo 
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Since the spectral radius r(Q\) is the only eigenvalue with positive eigenfunction for the strongly positive 
compact operator Q\ € ^(-SgJs ? ) by the Krein-Rutman theorem, we have r{Q\) = k(X). To satisfy 
d2.19t we assume that 

jfe(0) = / b{Q,a)e~ ai ^ D ^^ dr da > 1 (3.8) 
Jo 

and then choose Ao > such that fc(Ao) = 1. 

Summarizing what we have just shown and referring to Theorem |2.5| we can state: 
Proposition 3.1. Suppose d3.1| l- (l3.8l l. Then (Ao, 0) with fc(Ao) = 1 is a bifurcation point for the problem 

d a u + A(U(x),a)u + X/jl(U (x) , a)u = 0, a G J, a: 6 fl , 

u(0,x) = / b(U(x), a)u(a) da , x e , 
Jo 

Bu(a, x) = , a> 0, x & dfl , 

U(x) = I u(a, x) da , ie!!. 



There are Sq > and a unique branch {(A(e), u(s)) ; |e| < £o} of solutions emanating from (Ao, 0) with 
u(e) e L P (J, Wl B ) n Wp(J, L q ), u(e) ^Otfe^O, 

of the form 

u(e) = e( e- Ao J? Mo,r)dr e"^ Jo* op.')^ + z(e) ) ^ | £ | < £q (3 9) 

fiof/i A : (— £o,£o) — ► K + z : (— £o,£o) - ► L p (J, B ) (~l Wp(J,L q ) are continuous. Moreover, 
u(e)(a, x) > /or e S (0, £o) ant/ (a, i) G J x O. 

Remark 3.2. A /oca/ dependence of the data on u with respect to age is also possible. For example, one 
may apply Theorem 12. 5\ for diffusion terms of the form div x (D(u(a, x)) V u ) as well (see 1201 Ex.3.1] /or 
details). Moreover, the functions D and d in (13. It may also depend on x G Q provided the dependence is 
sufficiently smooth. In this case one needs Remark \4.2\ d) to verify (12.51 . 

3.2. ProofofTheorem ll.il Applying Proposition ^, ll to the problem 

<9 a w = div x (D(U(x),a)V x u) - Xp,{U(x),a)u , a E J , x <E n , (3.10) 

u(0,x)= f b(U(x),a)u(a)da, xEtt, (3.11) 

•/ o 

fot(a,a;) + (1 - S)d v u(a,x) = , a e J , x e dVL . (3.12) 

considered in the introduction, we obtain under the assumptions of Theorem 11.11 a branch of nontrivial 
solutions 

{(A(e), «(£)); |£| <£ } in R+ x (C*(J, C(fi)) n J, C(Q)) n C( J, C* 2 (ft))) , 

where the regularity of w(e) is due to standard parabolic regularity theory (e.g. see (3] Thm.9.2]), where 
j:=J\{0}. 

Let now S = in d3.12t and assume d 1.51 l, Then o\ = and, for any nonnegative solution (A, u) to 
( l3~T0l l-( [3~T2l , we have 

— / u(a, x) dx = — X [ fi(U(x),a)u(a,x)dx < ~Xfi(0,a) f u(a,x)dx, 
da in Jn Jn 

whence 

z{a) < z{Q)e- x t<>^ r)dr for z(a) := [ u{a,x) dx . 
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Moreover, by ( 13. lit and ( 11.51 ), 

z(0)= / / m 6(C/(a;),a)w(a,x)dada; < / " 6(0, a)z(a) da < z(0) / "' 6(0, s)e^^* r)dr da , 



n Jo 



and thus 

fc(A) = / 6(0, a)e~ X /o" "(0,r)dr dfl > j 



Since fc(A) is strictly decreasing in A and fc(Ao) = 1, we conclude A < Ao, that is, subcritical bifurcation 
occurs in d3.10t -( l3~T2l > in this case. This proves Theorem ll.il 

3.3. Subcritical bifurcation for Dirichlet boundary conditions. We consider an example involving Di- 
richlet boundary conditions. More precisely, let us consider 

d a u = D(U)A x u- \n{U{x),a)u , a e J , x e fl , (3.13) 

u(0,x)= f b(U(x),a)u(a)da , (3.14) 
Jo 

u(a,x) = 0, aeJ,xedn. (3.15) 

Note that the diffusion coefficients are independent of a G J. Suppose that D G C 1 (W^ B (f2), (0, oo)). If 
( 11.5b still holds, then we easily derive for any positive solution (A, u) of d3.13t - d3.15t analogously as above 
that 

z' (a) + <j\D(U)z{a) < — A^(0, a)z(a) for z(a) := / (p\{x)u{a, x) dx , 

Jo 

where (p\ is a positive eigenfunction to the principal eigenvalue u\ > of — A x subject to Dirichlet bound- 
ary conditions on dil. Thus, if in addition D(U) > D(Q) for < U G Wg B (Ct), then we deduce again 
that 

z(0)< I"™ b{<d,a)e~° lD{0)a e- x K » {0 ' r)dr daz{0) = k(X) z{0) . 



Hence A < Ao and thus subcritical bifurcation occurs also in this case. 

Observe that one may replace the diffusion term D(U)A x u in ( 13.131 l by div x (D(U (x))V x u) depending 
locally with respect to x on U and derive the same conclusion of subcritical bifurcation provided that 
oi (U) > cri(0)for0 < U G where <Jx{U) is the first eigenvalue of w h> -div x (D(U(x))V x w) 

and using a corresponding positive eigenfunction ipi — (pi (U) in the definition of z. 



3.4. An example with Holling-Tanner type nonlinearities. As noted in Remark 12761 one can also allow 
for ft, < in dl.61 >. We conclude with an example, which has been investigated in [ 14 1 in the case of linear 
diffusion: 

u 

d a u + A(U(x),a)u + u(U(x),a)u = Au± , a£j,xGQ,, (3.16) 

1 + u 



u(0,x)= b(U[x),a)u(a)da, x G Q , (3.17) 
Jo 

Bu(a,x) = 0, a>0,xedn, (3.18) 

U(x) = [ u{a,x)da, xeQ, (3.19) 



with A and B as in Subsection 13. II We impose the same assumptions d3.1t -( T3~7b as there, where we take 
q = p for simplicity. The strict positivity of /i(0, a) in d3.6l ) is not needed here. We also use the same spaces 
as in Subsection|3.1| 



Ei ■= Wis := Wl B (n) := {u G % 2 ; Bu = 0} «-». L p =: E , E, = %% 2/p 
Ei := L p (J, Wl B ) n Wp(J, L p ) and E := L p (J, L p ) . 
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Noticing that the "Holling-Tanner-type" nonlinearity can be written in the form 

u u 
1 + U _ U 1+u' 
problem d3.16t - (l3.19l l fits in the abstract form ( 11.6) . ( 11.7) by setting 

A(u,a):=A(U,a)+n(U,a)Tle£(W2 B ,L p ), h(a) := -1 , g( u ):= T ^—. 

Let V := (— 1, 1). Then ( 13.5b (with p = q), Remark l4~2l b) from the appendix, and [6, VII.Thm.6.4] ensure 
that the superposition operator of g (still denoted by g) belongs to C 1 (C( J, V p ), C( J, ~b^ P )), where we 
define V p := Wp~ B 2/p H (7(0, 7). Also note that g(0) = 0. Recalling that Ei ^ C(J, E q ), this implies 
g G C 1 (£, C(J, F)) for £ := B El (0, i?) with i? > sufficiently small and F := F f . To satisfy r(Q ) < 1 
(see Remark lZoT ). we assume that 

fc(0) < 1 (3.20) 

for 

fc(A) := f m 6(0, a) e (ATl)a e - # ^°^ dr e^ 1 K D(o,r)dr da ^ 
Jo 

where <j\ is the first eigenvalue of (— A x , B). Let ip\ G g be a corresponding positive eigenfunction. 
As in Subsection l3.1l we have Q\tp = k(X)(pi with 

Q A := f m 6(0,a)e( A T 1 ) a e -/ l >(0.Odr e / £ 'D(0,r)drA a;da; 

./o 

whence r(Q\) = fc(A), in particular, r(Qo) < 1 in view of (13.201 1. Thus we may invoke Theorem l2.5l and 
Remark [231 to conclude the existence of a branch of nontrivial solutions to (13.16) - d3~T9l ) emanating from 
the critical point (Aq, 0), where A > with fe(Ao) = 1. 

This local bifurcation result generalizes the (global) one of [ 14] in that nonlinear diffusion and nonlinear 
death and birth rates may be considered. However, we shall point out that in lfT4l a death rate depending 
on local position is considered (what can be considered in the present situation as well but requires some 
additional effort). 



4. Appendix 

We prove a result on the differentiability of superposition operators in Sobolev-Slobodeckii spaces that 
is used in the previous examples but might be of interest in other applications as well. The proof is similar 
to lfl8l Lem.2.7] or Prop. 15.4], where continuity properties are derived. 

To set the stage let £1 be an open and bounded subset of M", let V be an open neighborhood of in R k , 
and let / be a compact interval in R. Given a function f : I x V — > K define the superposition operator F 
off by 

F[u](a)(x) := f(a, u(x)) for u : Q, — > V and a G I , x e CI . 

We write / G C°' k ~ (I x V") provided that d% ~ 1 f 6 C(I x V) is Lipschitz continuous in x £ V uniformly 
with respect to a 6 J. 

Recall the definition of the norm in K fc ) (e.g. see Qj]): if q G (1, oo) and £ G (0, 1), then 

H u "W(n,R*) = l|u|l W*) + Ln l^i^" d(a! ' y) 

and if £ G (1,2), then 

71 
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Proposition 4.1. Let q G (n, oo), £ G (n/g, 2), anc/ 77 £ (0, £). Further, let V 6 := W|(f2, R fe ) n C(fi, V) 
£>e equipped with the Wfj -topology. Then F G C 2 ~ (1^,^(1, WJ(fi))) provided that f G C°' 4_ (7 x V). 
77ie Frechet derivative DF[u] at u G z'i g/ven fey 

(DF[«]fc)(a)(i) = d 2 f(a,u(x))h(x) , a€l,xeSl,he W^(n,R k ) . 

Proof. We may assume k = 1. 

(i) First, let £ G 1) and note that W^(Q) C(O). Fix neFj and choose an open neighborhood R 

of u(O) in V such that its closure R is compact and contained in V. Then, since / G C 0,3 ~(I x V), there 
is c (i?) > with 

\d%f{a,r)\<co(R) , rtR, ael, 



\d 2 f(a,r) - d 2 f(a, S )\ + \d 2 2 f(a,r) - d 2 2 f(a,s)\ < c (R)\r- 



r,s G R , a E I 



In the following we suppress the (fixed) variable a G I in / and its derivatives for the sake of readability 
and we set /' := d 2 f ■ Let h G W|(f2) with H^ll^/S^ sufficiently small so that u(f2) + fr(O) C -R. Then 

(F>]/i)(a)(z) := Q,/(a,«(x))ft(a;) = f'(u(x))h(x) 
by convention. The mean value theorem implies for x, y G Q: 

+ fe](a)(x) - FH(fl)W - (F>]/i)(a)(z) - [F[ u + ft]( a )( tf ) - F[u](a)(y) - (F'[u]/i)(a)(j/)] | 

[/'(«(!) + Th(x)) - f'(u(x))] dr (h(x) - %)) 



< 



1 r i 



JO 



/" (u(x) + arh(x)) da h(x) dr h(y) - / / /" (u(y) + ar%)) da dr h(y) 



Jo 



< c (R)\\h\\ oo \h(x)-h(y)\ + 



f"(u(x) + aTh(x)) dcrdr (h(x) - h(y)) h(y) 



Jo 



[f"(u(x) + cTTh(x)) - f"(u(y) + aTh(y))]dadrh(y)h(y) 
< c{R) \\hU \h(x) - h(y)\ + c (R) \\h\\l [\u(x) - u(y)\ + \h(x) - h(y)\] . 
Therefore, we obtain by definition of the norm in the Sobolev-Slobodeckii space W^(Cl) that 



\F[u + h](a) - F[u]{a) - (F'[u]h)(a) 



'wf (fi) 



<c(R)\\h\\% + c(R)\\h\\Uh\\l 



=(Wll2S{l 



Wq (H) 



Recalling the embedding W|(f2) ■— > C(f2) we thus deduce that 

\\F[u + h}- F[u] - F'[u}h\\ L ^ {Lwlm = 0(11^11^1^) , (h-tO), 

whence F : V$ — > Lao(I, W|(fi)) is Frechet differentiable at u G with derivative DF[u]h = F'[u]h 
for h G W£(Q). Moreover, since pointwise multiplication W|(fi) x W|(f2) — > W|(fi) is continuous due 
to 101 Thm.4.1] and £ > n/q, we have for v G W|(f2) with sufficiently small norm that (writing here and 
in the following e.g. f'(u) for the superposition operator at u induced by /') 



\DF[u + v] - DF[u] 



\c(W$(Cl),L aa (I,W$(ti))) - SU P 



||f (u + v)h - f (u)h\\ LaD{I w i m 



< c\\f'(u + v) - f'(u)\\ LaoXI W i m < c(R)\\v\\ w i {Q) 
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where the last inequality can be shown similarly as above (or also follows from [ 18, Lem.2.7]). This implies 

FeC !i -(V^L 00 (I,Wl{Sl))). 

(ii) Now let £ £ (1, 2) and rj G (1, £). Choose r G (n/q, 1) with r > rj — 1. Then pointwise multiplication 
W^{VL) x W| _1 (fi) ->■ V^' -1 (fi) is continuous, see g] Thm.4. 1]. Therefore, taking £ = r in (i) and using 
the chain rule we obtain for u 6 V£ and ft, G W|(fi) with ||^||^«(q) sufficiently small that 

\\F[u + h]- F[u] - F'[u]h\\ Loo{I>w?m 

< \\F[u+h}-F[u}-F'[u}h\\ Lao{I , Lq{n)) 



j2Uf(u+h)-r(u)-r(u)h)d : 



j u \\w v ^ 1 



3 = 1 



+ E II (/'(« + ^) "/'(«)) ^llwj-cn) 

^ °(ll^llw<(n)) + c l|/'( u + h ) ~ f'( u ) - f"i u ) h \\w-{n) H\wl(n) 

+ c||/'(m + /i) - ./"(u)llw-(o) ll^llwr^n) 
^ (IWIwftn)) + °(ll^llwj-(n)) iMLj(fi) + c \\ h \\w S {n) I^IUl(n) 
= °(\\ h Wwl(n)) > (^°) 

the last equality being due to the embedding W|(fi) ^ (fi). This shows that F : -> L^I, WJ(fi)) 
is Frechet differentiable at u e with derivative DF[m]/i = for h G W|(fi). 

Finally, to prove the Lipschitz continuity of DF choose £ G (77, £) and note that pointwise multiplication 
W| (fi) x W|(fi) ->• WJ(fi) is continuous due to g] Thm.4. 1]. Hence it follows for v G W| (fi) with 
sufficiently small norm that 

linen , 1 nr[ in \\f'(u + v)h - f{u)h\\ LUI ^ m 

\\DF[u + v]-DF[u]\\ c(w t (n)MIwSm) = sup r—— 

fcew|(n) 1 1 1 1 w| (o) 

< c ||/'(« + «) - f'(u)\\ LooiI ^ m < c(R) \\v\\ wHn) , 

where the last inequality stems from US Lem.2.7]. So F G C 2 ~ W^(fi))). The case £ = 1 is 

obvious. □ 

Remarks 4.2. (a) If£ G (n/q, 1) and / G C"^-(I x V), f/ien F G C 2 ~ (V^L^I, W|(fi))). 

Proof. See part (i) of the proof of Proposition ^. 11 □ 

(&)//£ G (1 + n/q, 2) and f G C°' 4 -(I x V), f/ien F G C 2 ~ (% £«,(/, W|(fi))). 

Proof. This follows exactly as in part (ii) of the proof of Proposition 14. ll bv observing that pointwise mul- 
tiplication Wg(Q) x W| _1 (fi) — > W| _1 (fi) for r G (n/q, 1) with r > £ — 1 is continuous according to 
H Thm.4. 1]. □ 

(c) For simplicity we refrain form taking into account an explicit dependence of f on x G fi, that is, we 
do not consider functions / :/ xflxF->Ein Proposition \4.1\ Such a dependence can be included 
provided that f (and its derivatives) are Holder continuous with respect to x, see [0 Prop. 15.4, Prop. 15.6]. 
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